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Abstract. The analysis of flare start-times confirms the periods found years ago 
(Aschenbach et al., 2004) in the near- infrared and X-ray light-curves related to the 
Sgr A* black hole. The assignment of the frequencies found to radial and vertical 
epicyclic frequencies v T , u v , respectively, as well as to the Kepler orbital frequency 
vk reveals resonances of v v /v T =7:2, and vk/v v =3:1. The highest observed fre- 
quency of 10 mHz is identified as the Kepler frequency corrected by the rotational 
frame-dragging frequency, as expected from the Lense-Thirring effect. These fre- 
quency assignments conclude a black hole mass of M = (4.10-4.34) xl0 6 M o and 
a spin of a = (0.99473-0.99561). 
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1. Introduction 2003) and X-ray (Aschenbach et al. 2004) 

light-curves of Sgr A*. 

The center of our Galaxy is suspected 

to contain a supermassive black hole. These candidate periods are suspect 

Observationally this conjecture is based on because of their enhanced power spectral 

the analysis of motions and orbit sizes of density in at least three Sgr A* obser- 

stars moving around some central mass vations, which include the October 26, 

(Schodel et al. 2002; Ghez et al. 2005; 2000 {Chandra, Baganoff et al. 2001), the 

Gillessen et al. 2009). These observations October 3, 2002 (XMM-Newton, Porquct 

suggest the presence of a fairly small vol- et al. 2004) and the near-infrared (NIR) 

ume containing a mass of several million observations of June 16, 2003 (Genzel et 

solar masses. Mass and spin of the black al. 2003). These observations were particu- 

hole have also been estimated from quasi- larly interesting because they showed fairly 

periodic oscillations (Abramowicz et al. large outbursts of Sgr A*. I have sum- 

2004; Aschenbach et al. 2004; Aschenbach marized the results, and I have proposed 

2004) which have been suggested to be candidate period-like structures, centered 

present in the near-infrared (Genzel et al. around 110 s, 219 s and 1173 s (Aschenbach 

et al. 2004). Later measurements, both in 
the NIR and the X-ray band, were not con- 
Send offprint requests to: B. Aschenbach elusive on these periods. In most of the 
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cases there was no indication of a period 
at all, but when a time structured signal 
was suggested by the data, it happened to 
be close to the candidate period of 1100 s. 
This includes a 1330 s period (Bclangcr ct 
al. 2006), and a -1500 s period (Meyer et 
al. 2008). However, each one of these ob- 
servations, taken as a single event, was not 
and cannot be considered to be of that sig- 
nificance which would justify a statistically 
unambiguous claim of the detection of a pe- 
riod. But the fact that these coincidences 
among the light-curve observations exist as 
far as the putative periods are concerned, 
is encouraging further study. 

The indication of more than one period 
in the light-curve data suggests the possi- 
bility that the light-curve is not dominated 
by just one period but that the signal is ac- 
tually modulated by one ore even more fre- 
quencies. So far, the analyses had to deal 
with light-curves which are usually subject 
to a large amount of noise (be it white, 
or red or even pink). Therefore I took the 
opportunity to analyse the starting times 
of a sequence of X-ray flares which were 
observed with XMM-Newton between 31 
March and 5 April 2007, which in principle 
provide an independent access to periodic 
patterns. 

2. Data and analysis 

In a recent paper Porquet et al. (2008) 
reported for the first time a high-level 
X-ray flaring activity of Sgr A* ob- 
served with XMM-Newton between 31 
March and 5 April 2007. Five flares in 
a row were detected. Their start-times 
were assessed in a quantitatively justi- 
fied way by Porquet et al. (2008). They 
were determined to 291913503, 292051530, 
292073530, 202084630, 292092330 in sec- 
onds relative to the on-board clock read- 
ings of XMM-Newton. The shortest separa- 
tion between any two of the flares is 7700 s, 
and 178827 s is the widest interval. Given 
these times I checked each interval between 
any two of the flares for accommodating, 
as closely as possible, a natural number of 



trains of a trial period covering the range 
from 10 s to 8000 s with a spacing of 1 ms, 
i.e. ~8x 10 6 periods were checked. The data 
provide 5 time-settings, or four indepen- 
dent time-intervals. The choice of four in- 
dependent time- intervals is, however, some- 
what arbitrary; one can use the time differ- 
ence between two consecutive events or one 
can choose the time difference measured 
against a reference, i.e. the first flare, the 
second flare, etc. I decided to take all pos- 
sibilities into account, which for five flares 
makes 10 measurements. It is obvious that 
these 10 time-interval measurements are 
not independent of each other, but this pro- 
cedure tends to reduce any biases in the 
measurements, e.g. large measurement er- 
rors for the one or other start-time. The 
problem with this approach is that this pro- 
cedure is likely to produce not only some 
solutions but also combinations of them, 
because of the oversampling of the infor- 
mation. But any eventual results can be 
screened for such events and eliminated af- 
terwards. 

The algorithm (sec Appendix 1) looks 
for a regular flare start-time pattern, keep- 
ing in mind that most of the flares have pos- 
sibly a brightness below the detection limit 
and are missed in the observations. The 
quantitative search determines the mini- 
mum residual between the observed inter- 
val and the time of cither n-timcs or (n+1)- 
times of the trial period. The values of the 
residuals are squared and added for the 
10 available intervals, divided by the num- 
ber of measurements (10) and after taking 
the square root they end up as root-mean- 
square (rms) value, which is then divided 
by the trial period. In other words, the 
mean relative rms-mismatch per measure- 
ment between observation and trial period 
is calculated. This is a dimcnsionlcss quan- 
tity, and this quantity can take values be- 
tween and 0.5. The inverse of that quan- 
tity I call 'goodness of fit' (GOF), which 
is the ratio of the trial period and the 
rms deviation from this trial period for the 
10 time-intervals. This is NOT a signifi- 
cance test, it is just a search algorithm for 
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Fig. 1. Pcriodogram of the trial periods computed from the X-ray flare start-times of the April 
2007 XMM-Newton observation. The ticmarks in the uppermost part of the graph delineate the 
positions or periods discussed in detail in the text. Data with a 'goodness of fit' value of less 
than 10 are not shown. For the definition of 'goodness of fit' see the text. 

Table 1. Periods in the X-ray flare start-times of the April 2007 XMM-Newton observa- 
tions. 



position 


period P 

(s) 


line width AP/P 
(FWHM) 


frequency v 
(mHz) 


mode assign. 


1 


100.016 


1.2xl0~ 4 


9.9984 


2^6 +V7 + f'J 


2 


102.777 


1.7xl0~ 4 


9.7298 


2l>6 + V7 


3 


109.980 


1.5xl0~ 4 


9.0926 


2V(, + IA) 


4 


188.046 


1.7xl0~ 4 


5.3178 


V% + "7 


5 


213.650 


2.6xl0~ 4 


4.6806 


V6 + "9 


6 


226.656 


2.6xl0~ 4 


4.4120 


El 


7 


1103.893 


1.4xl0~ 3 


0.90589 




8 


1569.126 


2.8xl0~ 3 


0.63730 


V 7 - Vg 


9 


3723.179 


7.5xl0~ 3 


0.26859 


Hi 



suspects! The 'goodness of fit' versus the trial period, which I call a periodogram, is 

shown in Fig. 1. 
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The periodogram shown in Fig. 1 re- 
veals nine prominent peaks which are listed 
in Table 1. The third column of Fig. 1 
shows the width of the outstanding peri- 
ods, which is very small, indicating that 
a search with a 1 ms stepping is too wide 
for any periods shorter than ^30 s. For the 
range with trial periods shorter than 35 s 
the screening was repeated with a step size 
of 10/xs. The corresponding periodogram is 
shown in Fig. 2. 

The peak at 33.3321 re-appears, and 
other peaks at shorter periods show-up. 
These large peaks, however, have in com- 
mon that the periods are exactly (to the 
fourth digit) given by Pi/n with n a natu- 
ral number. They do not provide indepen- 
dent information, but none of the periods 
with P>100 s, listed in Table 1, can be rep- 
resented as a simple linear relation between 
any one of the listed periods. 

However, when analysed in the fre- 
quency domain linear relations become ap- 
parent for P>100 s, and there are just three 
independent frequencies left, out of which 
the rest can be made up by linear combi- 
nations. Because the search algorithm is in 
the period domain rather than in the fre- 
quency domain these three frequencies pro- 
vide independent information. 

The nine periods consist of two groups 
with one group centered around ~106 s and 
a second group centered around ^214 s, 
and two periods at ^1104 s (v-;) and 
^3723 s (fg), respectively, as well as their 
difference frequency {v$). Both, the group 
of the three frequencies around ^106 s and 
the group of the three frequencies around 
~213 s each follow exactly the frequency 
relations for side band frequencies created 
by v-j and vg. In principle, there are eight 
side band frequencies for each group and its 
central frequency. Apparently, only three 
of these nine possible frequencies show up, 
for each group, with a high peak in the 
periodogram. The central frequencies v c 
of the two frequency groups are calcu- 
lated to ^ c .i=9.3611 mHz (P Cj i=106.825 s) 
and j/ Ci2 =4.68055 mHz (P c ,2=213.65 s). 
Obviously, ^ c ,2=^c,i/2, i.e., there is just one 



fundamental frequency, which is v c =is Ci i. 
In summary, the observationally fundamen- 
tal four flare time-intervals have provided 
three independent frequencies, i.e. v c (and 
v c /2), v 7 and v 9 . 

Of course, I have run a number of simu- 
lations assuming that five flares happen to 
occur at random over the entire observa- 
tion time of 230 ks. GOF values as high as 
observed are found, for one or another pe- 
riod, quite frequently for periods less than 
about 100 s, but they tend to occur much 
less frequent with increasing period, so that 
the periods close to 1000 s and above are 
outstanding. 

In addition, the periods showing up 
in the periodogram derived from the 
flare start-times are very much reminis- 
cent of the periods suggested earlier to 
be present in the Sgr A* light-curves be- 
cause of their enhanced power spectral den- 
sity, i.e. the observations of October 26, 
2000 (Chandra), October 3, 2002 (XMM- 
Newton) and the NIR flare on June 16, 
2003 (for a summary see Aschenbach et al., 
2004). The numbers to be compared are 
106 s (110 s), 214 s (219 s) and 1104 s 
(1173 s). The numbers without brackets 
are derived from the analysis of the flare 
start-times taken in 2007, and the number 
in brackets are from the light-curves taken 
before 2004. 

In view of this, I have not undertaken 
the tremendous task of calculating the 
probability that the observed periods are 
not random events, because that probabil- 
ity is definitely not zero. Any associated 
significance number just helps to decide 
whether to continue or stop the investiga- 
tions. Because of the astonishing overlap 
between the periods with those, which show 
some enhanced power spectral density in 
the flares mentioned above, it is useful to 
investigate the possibility that the periodic 
structures are associated with the Sgr A* 
black hole. 
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Fig. 2. Pcriodogram of the trial periods computed from the X-ray flare start-times of the April 
2007 XMM-Newton observation in the range of 10 s to 35 s using a screening of 10 /is steps. 
Data with a 'goodness of fit' value of less than 25 are not shown. 



3. Analysis and results 

The periodogram presented suggests the 
existence of three independent frequencies. 
Lacking any alternative model, which in- 
vokes more than just one frequency, I re- 
strict the following analysis to an attempt 
of assigning the observed frequencies to 
matter oscillations somewhere in an accre- 
tion disk which is bound by a rotating black 
hole. I suggest to consider that a hot spot is 
created from a pre-existing fluctuation by 
resonant waves traveling through the accre- 
tion disk. These waves are in the radial di- 
rection, transporting mass and energy from 
the outskirts to the black hole with a fre- 
quency v T , and in the vertical direction, 
mainly raising mass density and energy 
density, with a frequency v v . The amplifica- 
tion of a fluctuation is most effective if the 
waves get into resonance, i.e. v v / v T =n:m, 



with n and m natural numbers, and m 
as low as possible, preferentially m=l, or 
m=2. The resonance effect becomes even 
more dramatic if one or both of the waves 
get into resonance with the orbital mo- 
tion of the initial fluctuation, which is the 
Kepler frequency vk, i-e. VK./v v =\t:l, with 
k and 1 natural numbers, and 1 as small as 
possible, preferentially 1=1 or 1=2. Such a 
resonance can also be considered to hold for 
z/k and v T . The presence of parametric res- 
onances between the epicyclic frequencies 
has originally been suggested by Kluzniak 
and Abramowicz (2001a,b) to explain the 
3:2 high frequency quasi-periodic oscilla- 
tion pairs, which are observed in a num- 
ber of low mass X-ray binaries most likely 
containing a black hole companion. 

In the following I use the letters P for a 
period, u=l/P for the associated frequency 
and f2 for the frequency expressed in GR 
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units, i.e. c=G=M=l (c is the speed of 
light, G is the gravitational constant and 
M is the mass of the black hole) . v and fl 
are related through f2/V=27rGM/c 3 . 

In their 1973 paper Cunningham and 
Bardeen studied "The optical appearance 
of a star orbiting an extreme Kerr black 
hole. Recognizing that ".. the 'photon' tra- 
jectory can loop around the black hole 
any number of times', they concluded that 
there will be more than one image of the or- 
biting star seen by the infinitely remote ob- 
server, which are separated in time by one 
'photon' orbit, and they called these im- 
ages 'direct image', 'one-orbit image', 'two- 
orbit' image, etc. The effect on the ob- 
server is that the modulation of the light 
does not follow the Kepler orbital frequency 
but it is modified by the circular 'photon' 
orbit frequency f2 p h such that Oi/f2K=(l- 
^K/^ph)" 1 - is the frequency with which 
the light appears to be modulated in the 
observer's frame. Q p h is the same frequency 
as the rotational frame-dragging frequency 
(see Appendix 2). As long as OK/^ P h"Cl, 
ill w f2x, but if this condition is violated 
and the frequency ratio gets closer to one, 
there is a large difference and Q\ can be 
boosted up to several times of Qk- This is 
the case when the black hole has a spin 
of a ~ 1, and if the Kepler orbit radius 
is just a few gravitational radii or less. If 
this happens, the modulation of the light- 
curve is no longer with fix, and there is no 
frequency at the value of £!k- The by far 
highest frequency is then Q\. 

The frequency data shown in Table 
1 suggest that f r =fg, v y =vj and v\=v c . 
Using the relativistic expressions for Q p h, 
Qk, and Q r (see Appendix 2) a best- 
fit reveals Pk=1/^k=368 s if i/k/i/ v =3:1, 
and ^ v /^r=7:2. The value of P c was fixed 
at its measured value. The value of P v was 
fitted to 1085.3 s instead of the measured 
value of 1103.9 s, which is a difference of - 
1.7%, and P r was fitted to 3798.5 s instead 
of the measured value of 3723.2 s, which 
is a difference of +2.0%. The best-fit also 
provides the black hole mass and the black 
hole spin, which arc M = 4.22xl0 6 M Q 
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and a=0. 99519. The orbit resonance ra- 
dius is at r=1.49, which is 1.115 times 
the radius of the innermost stable circu- 
lar orbit or R = 0.93xl0 12 cm in physical 
units. An uncertainty of ±2% in the mea- 
sured values of v v and v x each, leads to an 
acceptable range for the black hole mass 
of 4.10xl0 6 M Q <M<4.34xl0 6 M Q and a 
range for the spin of 0.99472<a<0. 99561. 



4. Light-curves 

According to the frequencies found tempo- 
ral variations of the light-curves, both in 
the near-infrared as well in the X-ray band, 
should occur. Hot-spots are periodically 
generated at a rate of the vertical epicyclic 
frequency v w (every 1100 s). The occurence 
rate is, however, modulated by the radial 
epicyclic frequency v r =2i> v /7, such that 
the flare rate varies with i/ v x(l±f r /V v ). 
Accordingly, in snap-shot like observations 
periods can be found which vary between 
850 s and 1540 s among different observa- 
tions. 

Because of the resonance of vk and v v , a 
hot-spot, once created, has the chance to be 
re-heated periodically, and flare again. The 
time separation between such flares should 
be exactly equal to P v (1100 s), or 2xP v 
(2200 s). 

When the flares originate from two dif- 
ferent locations of hot-spots around the or- 
bit, the time separation between two con- 
secutive flares can differ from P v by as 
much as one orbital period, because of the 
optical transfer function of the black hole, 
which causes a delay between the hot-spot 
creation epoch and the time when the im- 
age passes the line of sight of the observer. 
The difference between the start-times of 
two consecutive flares is therefore some- 
where between 730 s<Atfl are <1470 s. 

At last, the light-curves will be modu- 
lated by the orbital motion of the hot-spot, 
but at a frequency of v\ (106 s) and ^ v /2 
(214 s). 
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5. Conclusion 

The periods indicated in earlier observa- 
tions of the near-infrared and X-ray light- 
curves from the accretion disk around the 
Sgr A* black hole are confirmed by the 
analysis of the flare start-times. The fre- 
quencies found imply resonances of 3:1 and 
7:2 suggesting that radial and vertical wave 
oscillations as well as the Kepler frequency 
arc involved. Taking into account the ro- 
tational frame-dragging, which should any- 
how be included in such an analysis (the 
Lense-Thirring effect, the Cunningham- 
Bardeen approach), the data are consis- 
tent with the presence of a supermas- 
sive, M = (4.10-4.34) xl0 6 M Q , and very 
rapidly spinning, a = (0.99473-0.99561), 
black hole. The inclusion of rotational 
frame-dragging, which is one of the key fea- 
tures of General Relativity, is mandatory in 
this analysis and, of course, interpretation. 

An open issue is the ad-hoc assump- 
tion that hot-spots are created from pre- 
existing fluctuations. The analysis pre- 
sented shows that the spin of the Sgr A* 
black hole is very close to a c —0. 0.9953, 
above which the anomalous orbit veloc- 
ity effect or 'Aschenbach' effect (Stuchhk 
et al. 2005) occurs (Aschenbach 2004). At 
a c =0.9953, the radius derivative of the or- 
bital velocity dv^ jdr = 0. This condition 
might cause an instability in the disk, like 
a shear-instability, which is being amplified 
by the radial and vertical waves. 
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Appendix 1 

A series of start-times of flares has been 
measured, and the time differences A; be- 
tween any two flares is searched for peri- 
odicities by scanning the values of Aj with 
a trial period P and a best-fitting natural 
number, which is cither n; or n; + l. The 
choice between ni and nj+1 is made by the 
lower residual. There arc N such intervals. 



a; = (A ; - m x P) 2 (1) 

6; = (A; - (ni + 1) x P) 2 (2) 

Cj = min(a;, b;) (3) 

x = v^EtiQ (4) 

The inverse of the root mean square 
value X divided by the trial period P is 
called 'goodness of' fit' (GOF). 

GOT = | (5) 



Appendix 2 

There are three cyclic modes associated 
with black hole accretion disks, which are 
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the Kepler frequency (J^k), the disk pertur- 
bation frequencies in vertical and radial di- 
rection called vertical (fly). Each frequency 
depends on the central mass M, the angu- 
lar momentum a and the radial distance r 
from the center. Equations [5] to [5] show the 
relations (Aliev & Galtsov, 1981). Equation 
[5] shows the dependence on the variables for 
the rotational frame-dragging frequency or 
f2 p h- The notation of c=G=l=M is used. 
Physical length scales are in units of GM/c 2 
and angular frequencies £1 are in units of 
c 3 /GM. r = 1 is defined as the gravitational 
radius r g . 

fl K = (r 3 / 2 + a)- 1 (6) 

n v = nl (i - ^ + (?) 
n% = nl (i - a + ^ - (8) 



! P h 2(l+Vl-a 2 ) ' 



(9) 



